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An important property of star-forming galaxies at z~l-2 is the high local star-formation intensities they maintain over tens of kilo- 
parsecs at levels that are only observed in the nearby Universe in the most powerful nuclear starbursts. To investigate how these high 
star-formation intensities affect the warm ionized medium, we present an analysis of the average spectra of about 50 such galaxies 
at za; 1.2-2.6 and of subsamples selected according to their local and global star-formation intensity. Stacking allows us to probe 
relatively weak lines like [Sn]/i/i6716,6731 and [Oi]/I6300, which are tracers of the conditions of the ISM and are undetectable in 
most individual targets. We find higher gas densities (hence pressures) in intensely star-forming regions compared to fainter diffuse 
gas and, overall, values that are comparable to starburst regions and the diffuse ISM in nearby galaxies. By modeling the Ha surface 
brightnesses and [Sii]/Hq- line ratios with the Cloudy photoionization code, we find that our galaxies continue trends observed in local 
galaxies, where gas pressures scale with star-formation intensity. We discuss these results in the context of models of self-regulated 
star foiTnation, where star formation determines the average thermal and turbulent pressure in the ISM, which in turn determines the 
rate at which stars can form, finding good agreement with our data. We also confirm the detection of broad, faint lines underlying Ha 
and [Nil], which have previously been considered evidence of either outflows or active galactic nuclei. Finding that the broad compo- 
nent is only significantly detected in stacks with the highest average local and global star-formation intensities strongly supports the 
outflow interpretation, and further emphasizes the importance of star-formation feedback and self-regulation in the early Universe. 
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1. Introduction 



The nature and evolution of galaxies is a result of the complex 
interplay between heating, cooling and dynamical processes in 
the interstellar medium (ISM). The cyclic, or competitive, na- 
ture of these processes determines the rate at which stars form 
and creates a feedback loop that determines the chemistry and 
structure of galaxies, along with their ISM. This interplay ulti- 
mately results in galaxies as we observe them today. 

In starburst galaxies, where the energy injection rate per unit 
volume from young stars is high, we may see the effects of the 
self-regulation of star formation, a process in which supernova 
feedback wi th its stron g energy injection may well play an im- 
ortant role ( Silkl2()0l]) - and drive the cloud vel ocity dispersion 
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.jJoung & Mac Lowll2006t 
bulent pressure (ISilkll2001l 



^asker & Brvan 2006) - as may tur- 
Blitz & Rosolowskv2006l) . At low- 



energy injection rates into the ISM, global shear is likely to play 
the most significant role in determining the peculiar velocities in 
dense massive clouds (Gammie et al. 1991). The observed tur- 
bulence is dominated by nonaxisymmetric gravitational instabil- 
ities at low star-formation intensity levels, like the few to several 
tens of km s"' seen in the Milky Way, but at high intensity lev- 
els mechanical energ y from the stellar population may play the 
most significant role (iJoung et al.i.2009: Ager tz et al...2009i') . 



* Data obtained as part of Programs ID 074.A-9011, 075.A- 
0318, 075.A-0466, 076.A-0464, 076.A-0527, 076.B-0259, 077.B-0079, 
077.B-0511, 078.A-0055, 078.A-0600, 079.A-0341, and 079.B-0430 at 
the ESO-VLT. 



Recently, we have proposed that the large Ha line widths 
observed in distant (z~l-3) intensely star-forming galaxies are 
driven by the mechanical energy liberated by young stars, a re- 
lationship which can be u nderstood within the context of self- 
regulated star formation (iLehnert et al.l 120091 : Le Tiran et al. 
201 1, in prep.). We hypothesized that mechanical energy is suf- 
ficient to keep the disk critically unstable against fragmentation 
and collapse, hence star-formation (Toomre Q~l). The intensity 
of the star formation in these distant galaxies is very high, similar 
to those in local starbursts but on a much larger physical scale, 
and it is maintained by large gas fractions and high mass-surface 
densities. These line wi dths do not appear to be driven by either 
cosmological accretion (iLe Tiran et al ]l2011h or gravitational in- 
stabilities. 

In this manuscript, we develop these arguments further 
through a stacking analysis of rest-frame optical emission lines 
of about 50 galaxies with redshifts of 1 .2 to 2.6 and high Ha sur- 
face brightnesses. Stacking allows us to analyze emission lines 
that are too faint to be observed in most individual galaxies 
such as [Sii]/l/l6716,6731 or [Oi]/16300, which are key indica- 
tors of the pressure and sources of ionization in the emission 
line gas. We analyze stacks of subsamples to investigate how the 
gas properties scale with star-formation intensity (rate per unit 
area). In a stacking analysis of a very simil ar sample, but focus- 
ing on trends with stellar mass and radius, IShapiro et al.l (|2009|) 
identified a broad line underlying Ho- and [Nii], but were un- 
able to differentiate whether this broad line emission was from 
outflows or active galactic nuclei. We confirm the detection of a 
broad component and argue that its presence in our stacks, where 
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the significance of the broad component increases with increas- 
ing star-formation intensity, can be interpreted as confirming the 
outflow hypothesis. 

2. Data analysis 

We use ESO archival data from a variety of programs with 
SINFONI on the ESO-VLT of a sample of more than fifty galax- 
ies in the redshift range 1.2-2.6; see Le Tiran et al. (2011) for 
further details. We have excluded all galaxies with recognizable 
AGN features in their data cubes - especially broad lines and 
high ratios of [NnJ/16583/HQ;. We produce for each object an in- 
tegrated spectrum using all the spectra in the data cube where Ha 
is detected above the 3-cr level, which are co-added after shift- 
ing Ha to its rest wavelength in order to remove any broadening 
due to the velocity field and weighted by the Ho' signal-to-noise 
ratio (SNR) of each object in order to maximize the final SNR. 
Uncertainties were measured using a Monte-Carlo method fit- 
ting 1000 realizations of the spectrum (Fig.[T]inset). 

To investigate whether the stacked emission line properties 
depend on the Ho- surface brightness (and/or redshift) of the 
galaxies used, we also made comparative stacks of objects in 
three bins: (1) z < 1.8 and star-formation rate (SFR) < 100 
Mgyr', (2) z > 1.8 and SFR < lOOMgyr' and (3) z > 1.8 and 
SFR > 100 Mpyr'^ All galaxies have similar isophotal sizes 
(iLehnert et alJ 120091 Le Tiran et al. 2011). To investigate the 
role that star-formation intensity might play in determining the 
characteristics of the integrated spectrum, for each bin and for 
the entire sample we also made one sub stack using only the 18 
brightest pixels in the Ha distribution of individual galaxies, and 
another using only the remaining, fainter pixels. Eighteen pix- 
els coiTespond to approximately the FWHM of the PSF (~0"6). 
Each brightest pixel stack contains about the same total fraction 
of the emission in its sample (~25%; Table[Tli. In Table[Tl Ngais 
is the number of galaxies included in the stack; <Esfr> is the 
average star-formation rate per unit area (Mq yr"' kpc"^), cal- 
culated using the c onversion factor fo r Ha luminosity to the star- 
formation rate from|Kennicut3(11998|); and Lna.brightest/LHa.whoie 
is the proportion of the Ha luminosity contained in the 18 bright- 
est pixels stacks compared to the stacks of whole galaxies. 

Table 1. Properties of the different stacked spectra 



bin 




<2sFR> 


LllQ'.brightest/LHQ'.whole 


all 


45 


0.6 


0.25 


(1) 


22 


0.2 


0.21 


(2) 


14 


0.6 


0.29 


(3) 


9 


1.0 


0.24 



The galaxies generafly have complex morphologies and the 
low spatial resolution of the data makes it difficult to accu- 
rately determine the relative centering and detailed morphol- 
ogy of the brightest pixels. They are usually within one reso- 
lution element but can have complex spatial distributions. High- 
resolution line imaging obtained using adaptive optics or HST 
continuum imaging show the morphologies of the highest sur- 
face brightness emission is frequently complex, cl umpy, and of- 
ten not at the dynamical or isophotal center (e.g. iLehnert et al.l 
l2009t iGenzel et al.ll201 iT) . Since all data sets reach roughly the 
same limiting observed Ha surface brightness level, over a wide 
range in redshift, we could not construct stacks at constant rest- 
frame surface brightness levels due to the strong impact of cos- 
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Fig. 2. [Sii]/l6716/[Sn]/l6731 line flux ratio as a function of elec- 
tron density (cm^-^), for a temperature of lO^K (blue line). Only 
results for the stacks of the brightest pixels of Ha emission are 
indicated (red crosses). The error bars indicate the Icr uncertain- 
ties in the measurements. For clarity, stacks with [Sii] line ratios 
below the low-density limit are not shown. The inset at the lower 
left shows the distribution of electr on densities n^ in the nuclea r 
region of nearby starburst galaxies dLehnert & Heckrnanll 19961) . 
as well as the mean electron density for the stack, including only 
the brightest pixels of all our objects (solid blue line) and the 
corresponding +lcr uncertainties (dashed blue lines). 



mological surface brightness dimming. The average rest-frame 
surface brightness of each stack is given in Table |2l 

3. Emission line properties 

We obtain log [Nii]^6583/Ha = -0.75 to -0.55 in all four stacks, 
and log [Sii]^^6716,6731 = -0.5 to -0.75. We did not de- 
tect [Oi]/I6300 in any of our stacks, with an upper limit of log 
[Oi]/Ha = -1.4 to -1.9 (Icr), similar to those in the integrated 
spect ra of nearby star-forming galaxies (Lehnert & Heckm^ 
11994 and Hn regions. The flux ratio of [Sii]/I6716/[Sii]^6731 
is about 1.2 to 1.4, near the low-density limit (1.45; Fig. |2]i, 
and the highest values are found in the stacks with the highest 
surface-brightness Ha emission. For most of the stacks the den- 
sity must be very low (n^ ~ 10-100 cm"^, see Fig. [2]), whereas in 
the regions of the highest Ha surface brightnesses in each galaxy 
we find values of about 100 to 500 cm"-', with a mean of ~200 
cm"^ (Table |2j Fig.|2]l. These values are similar to nearby star- 
bursts, which show strong evidence of driving energetic outflows 
(Fig. H ILehnert & Heckmanll 19961) . Finding that the full stacks 
and the stacks excluding the brightest pixels are at (or near) the 
low-density limit suggests that the most extended emission con- 
tributes most of the flux and indicates that these outer regions are 
more like the diffuse interstellar medium in ne arby disk galaxies 
(ILehnert & Heckmanll 1994t I Wang et aOl998h . 

We also find no significant differences in the widths of the 
various lines analyzed. Although there is a slight tendency for 
[Nn]/i6583 to be systematically broader than Ha in all stacks, 
all individual width measurements are the same within the un- 
certainties. We do find a trend for the narrow components of 
both lines to be broader for the stacks with high average surface 
brightnesses. This is related to the trend for the most intense star- 
forming regions to have the broadest lines (Lehnert et al. 2009]). 

Several of the spectral stacks have an apparent broad line 
component underlying the region around Ha (Fig. [3]), as noted 
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Fig. 3. Zoom of the region of Ha and [Nn] for 12 different stacks 
(black lines). From left to right: the integrated spectrum of all the 
emission from each galaxy (whole), only the brightest 18 pixels 
in Ho- (brightest), and all other pixels (faint); from top to bottom: 
all spectral bins, bins (1), (2), and (3). The dashed blue line is the 
best fit to the strongest component of each line, the green dashed 
line the fit to the weaker offset component, and the red line is 
the total of the fitted lines (see text for details). For each stack 
we show the residuals to the fit underneath each spectrum (in 
blue for stacks that do not require a broad component and in red 
for those that do). Velocity offsets range from ~100-500 km s"' 
with the highest value seen in the stack with the highest average 
surface brightness. 



by IShapiro et al.l (l2009t) . It is not our intent to accurately con- 
strain the properties of this component but to robustly demon- 
strate that some of the data are better fitted by including a broad 
line. To this purpose we adopted a Bayesian approach, as analy- 
sis of such data based on minimization is often not very con- 
straining due to the difficulty of accurately estimating uncer- 
tainties (see Table |2]i. We computed the relative strength of ev- 
idence, Bayes factor B, for a model that includes a broad line 
versus one without. Generally speaking, ln{B) > 5 is considered 
to be strong evidence in favor of th e broad-Une models (inter- 
preted against the Jeffreys' scale; see lTrottal2008l and references 
therein). We find strong evidence in favor of broad lines in the 
stacks with the highest average surface brightnesses (see B val- 
ues listed in Tab.|2]i. In general, our fits (a s derived from model s 
with the highest values of B) agree with IShapiro et al] (l2009h . 
The stacks of the brightest 18 pixels have the lowest signal-to- 
noise as each is composed of only ~25% of the total flux and yet 
it is those stacks where the need for an additional broad compo- 
nent is the most robust. We also find that the best fits are broad 
(FWHM~1500 kms"') and that the total flux ratio of the broad 
to narrow Ha emission lines is about 30% for those stacks with 
a significant broad feature. 

4. Discussion and conclusions 

The results of the broad line analysis (§ O suggest a close 
relationship between the source of the broad emission and 
the star-formation intensity in these galaxies. Nearby galaxies 
whose starbursts are as intense as those observed in our dis- 
tant sample show strong evidence of driving large-scale out- 



flows (iLehnert & Hecknianlll996t) . In fact, one can find close 
analogs between our spectra with or without a broad component 
and th e low-redshift starburst galaxies in Lehnert & Heckma^ 
(Il995h . sugges ting both samples mu st have similar phenomenol- 
ogy. Although IShapiro et al.l (|2009|) have already suggested that 
broad lines in stacked spectra of z~2 intensely star-forming 
galaxies may indicate outflows, since they stacked according 
to stellar mass and radius (comparing nuclear and off-nuclear 
stacks) they were unable to rule out that the broad lines were due 
to AGN. The trends we see with average star-formation intensity 
contradicts the AGN hypothesis, simply because many intensely 
star-forming regions in these galaxies are off-nucleus, whereas 
the nuclear regions often hav e very low surface brightness (e.g. 
iForster Schreiber et al.ll2009l) . 

In the outflow scenario we expect, in addition to single 
narrow components of Ha and [Nn]/l/l6548,6583, three more 
lines of Ha and [NiiJ that are offset and likely very broad (e.g. 
iLehnert & HeckmanI Il996h . By making five physically moti- 
vated assumptions we avoid introducing nine additional free pa- 
rameters: (1) all three additional lines have the same offset veloc- 
ity, (2) velocity dispersions are the same for each Ha and [Nii] 
line component and are equal to the offset velocity, (3) the flux 
ratio of [Nii]/Hq' is given by fast shock models (ratios ranging 
fro m ~0.2-l for the velocities we considered in this modeling, 
see lAllen et al.l20d8h . (4) the velocity offset and the shock speed 
are the same, and (5) the flux of the offset Ha component is 10% 
of that of the main Ha line. This last value provides a match to 
the estimated peak fluxes in the best-fit single broad component. 
This reduces the number of free parameters to only one, the off- 
set velocity, while the other eight are constrained. 

We find that these fits are as significant as assuming a sin- 
gle broad component for velocity offsets of a few 100 km s 'and 
narrow-to-broad Ha flux ratios of 10%. These values are similar 
to those in the extended (i.e ., wind) emission in nearby starbursts 
(ILehnert & HeckmanI 1996h . The derived pressures are also sim- 
ilar to those in nearby starbursts, which is a necessary condition 
for driving winds (Fig.|2]i. There is also enough mechanical en- 
ergy in the shocks to power these flows. The mechanical energy 
output at a SFR of 150 Mq yr"' (typical of our high surface 
brightness sample) is about 10^"* erg s"'. The fraction of the to- 
tal energy fr om shocks fast en ough to explain the broad lines is 
about 1-2% jAUen et al.ll2008l) . If we assume that the mechani- 
cal energy output from stars is efficiently thermalized, we only 
require around 1-10% of the mechanical energy to energize the 
broad, blueshifted line emission. 

What do these results imply about the nature of the ISM in 
these high-redshift galaxies? We have already argued that the 
galaxies have high pressure. Figure |4] illustrates that the warm 
ionized medium in nearby star-forming and starburst galaxies 
and in our galaxies forms a continuity - a one-parameter fam- 
ily. Going from low to high Ha surface brightness we progress 
from diffuse ISM in nearby galaxies, through Hii regions (and 
their surroundings) to the nuclei of nearby starburst galaxies 
(e.g. lWang et"anil998l) . On average our galaxies lie at the high 
surface-brightness, low [Sii]/Ha end of the relationship, similar 
to the positions of local powerful nuclear starbursts. 

This continuity can be explained through a range of ISM 
pressures and radiation field intensities. We mo deled the data 
as ph otoionized clouds (using the code Cloudy; iFerland et al.l 
using the ionizing spectrum of a young (10^ yrs) stellar 
population forming stars at a constant rate with a Salpeter IMF 
(Leitherer et al. 1999j) and a range of column and volume densi- 
ties. Our results suggest that the galaxies have high average gas 
densities (~10-few x 100 cm"^) and high ionization parameters 
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log Hjj (cm '}= 3 
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<i;,i„> (erg s"' kpc"^) 

Fig. 4. [Sii]/l/l 6716,673 l/Ha flux ratio versus average Ha sur- 
face brightness for our 12 different stacks (see Tab.|2| of galaxies 
(red circles, with 1-cr uncertainties) superimposed on data for the 
diffuse emission in nearby star-forming and starburst galaxies 
from Wang et al. ( 1998) ( black dots). The line s represent Cloudy 
photoionization models (iFerland et al.lll998h with various col- 
umn densities (10^", 10"' and 10"-^ cm in green, red and pur- 
ple, respectively), ionization parameters (U=10"^ to 1 from left 
to right along each line) and densities (1,10,100,1000 cm"-' as 
solid, dashed, dot-dash and dotted, respectively). The legend at 
the upper right indicates the colors of the lines for different col- 
umn densities, while the legend at the lower left indicates the 
line styles for each volume density. 



(log U~ -2 to 0, where U is proportional to the relative inten- 
sity of the photon field di vided by the total gas density). This 
supports the hypothesis of IWang et al.l (Il998h that this relation- 
ship can be understood as an underlying proportionality between 
the thermal pressure and the mean star-formation intensity in a 
photoionized gas. 

What is the source of this underlying proportionality? Star 
formation might be regulated by the average pressure in the ISM 
(ISilklll997h . which IWang et all ( Il998h suggested may itself be 
either regulated by the mechanical energy injection from star 
formation or related to the hydrostatic or turbulent pressure. 

If the mechanical energy from massive stars is controlling 
the over-pressure, we would expect the pressure to increase lin- 
early with the star-formation intensity. This can be estimated as 
Pgas=M''~E''^ R*~ (where Pgas is the gas pressure, M the mass 
loss and entrainment rate, E the mechanical energy output, and 

the radius over which the en ergy a nd mass output occurs, 
e.g. lStrickland & Heckmanll2009h . From lLeiflierer et al.l (Il999l) . 
we can estimate the mechanical energy and mass output rate 
from star formation. Adopting an equilibrium mass and energy 
output rate for continuous star-formation over 10^ yrs, we esti- 
mate pressures of 1x10""^ dyne cm"^ (one-sided) for 0.5 Mq 
yr"' kpc"^. If the mass entrainment rate is a factor of a few, 
this predicted pressure would be some what higher (consistent 
with that observed in M82 for example; I Strickland & HeckmanI 
I2OO9I) . For our stacks, the results of the photoionization models 
suggest thermal pressures of 6x10""^ to 6x10"" dyne cm"^. 

In an alternative interpretation (which may be particularly 
appropriate at low star-formation intensities), the pressures are 
set by gravitational processes. Star-formation intensity is related 



3/'^ 1/2 

to a cloud-cloud collis ion model, which gives Ssfr ^gas^totai 
(ISilk & Normanll2009h . The pressure in the ISM can be related 
to gravity or turbulence through Pgas - PgasO-gas = ^GEgasStotai, 
where Pgas and pgas are the gas pressure and density, respectively, 
and G is the gravitational constant. Combining them suggests 
that ZsFR P(2gas/5^totai)'^^- We can normalize this relationship 
empirically with the appropriate values for the Milky Way - ISM 
pressure of ~3000 K cm"-', star-formation intensity ~2.5xl0"-' 
Mq yr~' kpc"^ and a gas fraction of about 10%. The gas fraction 
in distant galaxies is likely to be higher, a few 10s of percent 
(iDaddi et al.ll2010l) . Using the average pressure from our pho- 
toionization modeling and the relation between star-formation 
intensity and pressure based on the Milky Way scaling suggests 
that SsFR ~1 M0 yr"' kpc"^. 

If star formation is driving the gas pressure, then it is 
likely that turbulent pressure is similar or even higher than 
thermal pressure, especially at high star-formation intensities 
(iJoung et al.ll2009l) . In our stacking analysis, we find typical Ha 
line velocity dispersions of about 125 km s"' and densities in the 
warm ionized gas of about 30-300 cm"-'. If the clouds formed 
through turbulent compression, the pre-shocked material would 
have turbulent pressures; Ptui-b=PgasO"eas~ 10"'' to 10"'° dyne 
cm"^. This is higher than the thermal pressures we observe in the 
emission line clouds of nearby galaxies, as suggested by models 
of ISM pressure r egulated by feedback from star formation (e.g. 
iJoung et al1l2009l) . but in the range of what we find at z~l-2. 

In summary we find that these distant galaxies have high 
ISM pressures and drive outflows, at least at the highest average 
Ha surface brightnesses. These are two characteristics similar 
to intense starbursts at low redshift. Overall, we favor a picture 
where the pressure in the ISM is determined by the intensity of 
the star formation and where feedback sets the scaling between 
them. Since the pressure is being regulated and it likely deter- 
mines the nature of star formation, this suggests that the star 
formation is self-regulating in these high-redshift galaxies (ISiUd 
[T997I) . This extends our earlier conclusion that star formation 
is likely to regulate the pressure of the ISM based on spatially 
resolved properties of individual distant galaxies with high Ha 
surface brightness (Lehnert et al. 2009; Le Tiran et al. 201 1). 
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Fig. 1. Star-formation rate as a function of redshift for all galax- 
ies used. Red dashed lines correspond to the boundaries between 
the 3 different bins used to make comparative stacks, which are 
labeled (1), (2), and (3). As an inset on the left, we show the 
stacked spectrum (signal-to-noise of Ho' weighted average) of all 
the galaxies in the sample. The lines of He, [Nii]/l/l6548,6583, 
and [Sii]/(/l6716,6731 are all significantly detected, but we can 
only set an upper limit on the [Oi]/I6300 emission. 
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21.3 


1.3±0.2 


1 77+107 
' ' '-86 


all-faintcr 


25 9±0 5 


40.8 


244±3 


242 


5.8±1.7 


2030±870 


-270±340 


0.2 


11.5 


1 7±0 2 






23.4±0.6 


40.3 


241±5 


212 










-0.1 


1.6±0.3 






5 6±0 2 


40.6 


270±8 


230 










-0.1 


1 2±0 3 


981+--^' 




I o.o±U.O 


40.3 


234±5 


207 










-1.6 


1 7±0 4 




24 7±0 8 


40.9 


236±5 


231 










5.0 


1 7±0 3 




(2)-brightest 


7.2±0.3 


41.2 


265±8 


247 










4.6 


1.2±0.3 


184+793 


(2)-fainter 


18.0±0.7 


40.8 


227±6 


224 










2.4 


1.9±0.5 


(3) 


57.0±1.3 


41.1 


264±4 


278 


19.9±3.2 


1400±290 


-150±1I0 


0.3 


22.6 


1.6±0.3 




(3)-brightest 


13.5±0.4 


41.5 


292±7 


322 


7.1±0.8 


1160±160 


15±68 


0.5 


26.8 


1.2±0.3 


190+201 

' -141 


(3)-fainter 


44.7±1.2 


41.1 


259±5 


264 


12.5±3.4 


1550±630 


-270±220 


0.3 


8.5 


1.7±0.4 





Table 2. Detailed properties of all the stacks. Column (1) - The four bins used in the analysis are all objects, and bins (1), (2) and (3), see §[2] 
Tab.[T] and Fig.[T]for details. Within each bin, -brightest refers to stacks made using only the brightest 18-pixels area (which roughly corresponds 
to the area of one seeing disk) of Ho- emission per object, while -fainter refers to stacks made using only the remaining, less intense Ha emission 
pixels per object. Column (2) - Ha luminosity of the stack in units of 10*' erg s"'; Column (3) - The logarithm of the average Ha surface 
brightnesses of each stack in units of erg s"' kpc"^. <EHff> is calculated by dividing the total Ha luminosity of the stack by the total isophotal 
area of the galaxies making up the stack. The total isophotal area is estimated including all pixels that are 3 times the RMS of the background 
of each data set (see Lehnert et al. 2009 and Le Tiran et al. 2011 for details). Column (4) - FWHM of the Ho- line obtained without weighting 
the widths of individual objects by their Ho- flux (in km s"'); Column (5) - Ha FWHM obtained when weighting the widths of individual objects 
by their Ha flux (in km s"'); In the next columns, 5-8, we provide the characteristics of our fits of a single broad component underlying the 
Ha and [Nn] emis sion. However, we emphasize that we do not think of these fits as realistic (see text for details) but are given for comparison 
with lShapiro et al.l 1(20093. The parameters are only provided for fits that have high significance. Note the considerable uncertainties in the fits, 
suggesting they are not very constraining, which is the reason we adopted a Bayesian approach; Column (6) - Ha luminosity of the single broad 
line component fitted to the spectrum (in units of 10*' erg s"'); Column (7) - FWHM of the single broad component fitted to the spectrum (in km 
s"'); Column (8) - velocity offset between the best-fit broad component and the narrow Ha line component (in km s"'), where negative numbers 
indicate a blue-shifted broad component; Column (9) - flux ratio of the broad to narrow Ha line components; Column (10) - Bayesian likelihood, 
ln(B), of the model with a broad line compared to the model without, where values above 5 indicate a significant enhancement of the fit quality 
by including a broad line and negative values rule out a broad component; Column (11) - [Sii]/i6716/[Sii]/i6731 line flux ratios; Column (12) - 
electron densities (in cm"') derived from the [Sn] ratios given in Column (11), with their Icr uncertainties. Throughout our analysis, we adopt the 
cosmology Ho =70 kms"' Mpc"', fl^ = 0.7 and = 0.3. 



